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We present an analysis of the quantum adiabatic algorithm for solving hard instances of 3-SAT 
(an NP-complete problem) in terms of Random Matrix Theory (RMT). We determine the global 
regularity of the spectral fluctuations of the instantaneous Hamiltonians encountered during the 
interpolation between the starting Hamiltonians and the ones whose ground states encode the so- 
lutions to the computational problems of interest. At each interpolation point, we quantify the 
degree of regularity of the average spectral distribution via its Brody parameter, a measure that 
distinguishes regular (i.e., Poissonian) from chaotic (i.e., Wigner-type) distributions of normalized 
nearest-neighbor spacings. We find that for hard problem instances, i.e., those having a critical ratio 
of clauses to variables, the spectral fluctuations typically become irregular across a contiguous region 
of the interpolation parameter, while the spectrum is regular for easy instances. Within the hard 
region, RMT may be applied to obtain a mathematical model of the probability of avoided level 
crossings and concomitant failure rate of the adiabatic algorithm due to non-adiabatic Landau- Zener 
type transitions. Our model predicts that if the interpolation is performed at a uniform rate, the 
average failure rate of the quantum adiabatic algorithm, when averaged over hard problem instances, 
scales exponentially with increasing problem size. 

PACS numbers: 



I. INTRODUCTION 

Can quantum computers solve NP-complete prob lems 
using physical resources of time, space, and energy that 
are all bounded by polynomials in the size of the prob- 
lem? Most computer scientists are skeptical of such a 
possibility 0. However, recently a more physics-inspired 
perspective has arisen that is causing some to rethink 
this question. In 2001, Farhi et al. presented a quantum 
adiabatic algorithm for solving an NP-complete problem, 
and showed via numerical simulations on a sequence of 
progressively larger problem instances that the running 
time of this algorithm appears to grow only as a polyno- 
mial in problem size 0- By contrast, all known classi- 
cal algorithms for solving NP-complete problems require 
a running time that scales exponentially with problem 
size in the worst case 0- If the polynomial scaling of 
the adiabatic algorithm is correct, this would represent a 
monumental result for the field of quantum computing, 
as it would bring a host of useful but hard computations 
within the domain of computational tasks that can be 
performed exponentially faster on quantum computers 
than classical ones. With such extraordinary promise, it 
behooves us to understand the adiabatic algorithm in full 
detail. Unfortunately, it has proven to be exceedingly dif- 
ficult to obtain analytic results on the scaling behavior of 
the quantum adiabatic algorithm. Instead, for the most 
part, researchers have relied upon numerical simulations 
of small problem instances, typically involving 23 vari- 
ables or less 0- By comparison, a modern-day classical 



algorithm for solving 3-SAT can routinely solve instances 
containing several thousand variables 0, • It is ques- 
tionable whether numerical results based on 23-variablc 
simulations can be extrapolated reliably to 2000- variable 
instances. Thus a more analytic approach is needed. 

In this paper we develop such an analysis based on 
Random Matrix Theory (RMT) HUH, which is a sta- 
tistical description of complex quantum systems in which 
detailed knowledge about particle interactions is aban- 
doned in favor of a description in terms of random in- 
teractions. Such a description is typically found to be 
applicable to complex Hamiltonians without fundamen- 
tal symmetries, and with very large phase space (high 
dimension). In this case, the level density p(E) and 
the nearest-neighbor spacings (NNS) distribution assume 
universal laws, and only the statistical properties of the 
levels are of interest. Such Hamiltonian systems are usu- 
ally highly irregular, disordered, and chaotic. 

RMT has proven to be a successful method for pre- 
dicting properties of complex quantum systems that look 
superficially very different in terms of their energy eigen- 
spectra. By working with the nearest-neighbor level spac- 
ing fluctuations rather than the raw eigenspectra, deep 
similarities between apparently different physical systems 
have been revealed, and several hard to calculate prop- 
erties, such as transition rates, have been determined. 
In essence, a tractable model Hamiltonian can be used 
to make predictions about an intractable one provided 
their energy spectra can be described by similar NNS 
distributions. 
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In standard RMT, Hamiltonians are drawn from an 
ensemble of orthogonal or unitary Gaussian matrices 
(the Gaussian orthogonal or unitary ensembles, GOE 
or GUE), which corresponds to Hamiltonians with in- 
teractions of all possible particle ranks d (e?-body inter- 
actions). In other words, every element of Hilbert space 
is assumed to be connected to every other by an inter- 
action strength that is given by a random number. The 
level density corresponding to such a model is the cele- 
brated "semi-circle law" of Wigner Q , which predicts a 
power dependence between the density of states and en- 
ergy Physical complex quantum systems (such as large 
atoms or nuclei) on the other hand, are better described 
by interactions of rank two, that is, two-body interac- 
tions at the most. Such Hamiltonians can be obtained 
from the so-called "two-body random ensemble" (or 
embedded GOE) that well describes complex nuclear and 
atomic level spectra, but gives rise to a Gaussian, rather 
than polynomial, density of states. For Hamiltonians 
drawn from either a GOE or an embedded GOE, the 
distribution of nearest-neighbor level separations is uni- 
modal with a long tail, known as the Wigner distribu- 
tion 0, 0]. Such distributions are typical for complex 
systems without symmetries, which results in highly ir- 
regular energy spectra and chaotic dynamics. Conversely, 
Hamiltonians corresponding to physical systems subject 
to symmetries and conservation laws typically display 
regular energy spectra and Poissonian nearest-neighbor 
level separations. Such distributions fall off much faster 
than the Wigner distribution, and decay monotonically. 
The Brody distribution interpolates between those two 
distributions with a single Brody parameter q, where the 
limit g = corresponds to the Poissonian limit while 
g=l gives the Wigner distribution. 

A priori, the Hamiltonians arising in computational 
problems such as fc-SAT appear to have such a special 
structure that they are unlikely to be described by ran- 
dom interaction matrices. We can assess this by char- 
acterizing the irregularity of the NNS distribution of the 
instantaneous Hamiltonian of the adiabatic algorithm. If 
we find a Brody parameter close to zero, RMT cannot 
be used, while a Brody parameter closer to one indicates 
that RMT theory can predict global properties of the 
Hamiltonian dynamics reliably. Note that the scaling of 
the level density itself (polynomial for GOE, exponential 
for embedded GOE) is irrelevant for this determination. 

We determine the NNS distribution of the instanta- 
neous Hamiltonians solving 3-SAT problems by generat- 
ing random soluble problem instances (with exactly one 
solution) with a fixed ratio of clauses to variables, deter- 
mining for each the eigenvalue distribution, from which 
the fluctuations can be obtained. In short, the results re- 
veal a systematic change in the spectral regularity of the 
instantaneous Hamiltonians during the course of the adi- 
abatic algorithm. In the initial phase of the interpolation 
for especially hard problem instances, the statistical NNS 
fluctuations conform to a regular, Poisson-type distribu- 
tion. Later in the interpolation, the fluctuations conform 




FIG. 1: Avoided level crossings can occur when adjacent lev- 
els assume the geometry of converging hyperbolae. Both the 
slope difference Am and minimum gap AE affect the proba- 
bility of a transition occurring. 

to an irregular Wigner-type distribution instead. We also 
find that irregular spectra only occur for computationally 
hard problem instances. 

In this paper, we predict the scaling of the failure rate 
of the adiabatic algorithm for a fixed ratio of clauses to 
variables at a given point in the interpolation process, as 
larger and larger problem instances are considered. The 
adiabatic algorithm fails when the system spontaneously 
transitions from its ground state into any excited state. 
If we make the conservative assumption that the only 
source of non-adiabatic transitions are of the Landau- 
Zener (LZ) type [ill, i.e., localized transitions between 
adjacent levels at avoided crossings where the energy lev- 
els locally assume the geometry of hyperbolae, then we 
can obtain a lower bound on the transition probability 
using RAIT. Additional failure modes can only make the 
failure rate of the algorithm worse. 

With these assumptions and the use of RMT, we can 
then determine the transition rate from the ground state 
averaged over an ensemble of problem instances having 
the same ratio of clauses to variables. There are two 
model-dependent quantities in this result: the average 
ground state level spacing and the typical size of LZ 
asymptotic slopes. Ultimately, these quantities are re- 
lated to the parameters characterizing problem instances 
of the type being solved, i.e., the ratio of clauses to vari- 
ables. Hence the transition rate, at a given point in the 
interpolation, is related to the difficulty of the problem 
instances at that point. The rest of the paper is organized 
as follows. Section II describes the quantum adiabatic 
algorithm and the Landau-Zener transition probability. 
Section III summarizes prior research on the scaling prop- 
erties of the algorithm, i.e., whether the adiabatic cri- 
terion can be met if the interpolation is completed in 
polynomial time. Section IV introduces the concepts of 
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Random Matrix Theory and Landau-Zener transitions 
needed for our analysis. Section V reports on our nu- 
merical experiments and the gap fluctuation phenomena 
they reveal. Section VI uses the phenomena to justify a 
random matrix analysis of the adiabatic algorithm and 
describes the implications on the scaling of the quantum 
adiabatic algorithm. We discuss the distribution of gap 
energies in an Appendix. 
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where e«l. If this criterion is met, we can be sure the 
system will evolve into the desired state. But it is not im- 
mediately clear how quickly we can interpolate between 
and while ensuring this adiabaticity criterion is not vio- 
lated. 



II. THE ADIABATIC ALGORITHM 



The idea behind the quantum adiabatic algorithm is as 
follows. If a quantum system is prepared in the ground 
state of a time-independent Hamiltonian Hq, and if wc 
then cause the Hamiltonian to change from Ho to a final 
form Hi in T steps, e.g., by driving it linearly 

H{^) = (l-|)F + |ffi, (1) 

then the Adiabatic Theorem of Quantum Mechanics 01 
guarantees that the system will remain in the ground 
state of the instantaneous Harniltonians H(t), provided 
the change is made sufficiently slowly, i.e., adiabatically. 
Thus, if the final Hamiltonian can be made to encode 
a computational problem such that the ground state of 
corresponds to the solution to this problem, then the nat- 
ural quantum mechanical evolution of the system under 
the slowly changing Hamiltonian H(t) would carry our 
initial state into a final state corresponding to the solu- 
tion. A final state measurement would then reveal the 
solution. The key question is how quickly can one drive 
the interpolation between the initial and final Harnilto- 
nians while keeping the system in the ground state of 
the instantaneous Harniltonians passed through. If the 
shortest feasible interpolation time scales polynomially 
with increasing problem size, the quantum adiabatic al- 
gorithm would be deemed "efficient" , otherwise it would 
be deemed "inefficient" . An alternative way of stat- 
ing this is to ask under what conditions the passage from 
H and Hi can be performed adiabatically 0. If the 
minimum eigenvalue gap between the ground state Eq 
and first excited state Ei of the instantaneous Harnilto- 
nians is given by <7 m i n , where 

g m m=mm T [Ei(t)-E (t)} (2) 

and the matrix clement between the corresponding pair 
of eigenstates is 

(^)i,o = {E 1 ;t\^\E ;t), (3) 

then the Adiabatic Theorem asserts that the final state 
will be very close to the ground state of Hi{T) , i.e., 

|(i?o;ThA(T))| 2 > l-e 2 (4) 



III. PRIOR ANALYTIC RESULTS 

Prior analytic studies of the adiabatic quantum algo- 
rithm have yielded mixed results. Farhi ct al. analyze 
several models in which the gap behavior can be com- 
puted analytically, and be shown to decrease polynomi- 
ally in the problem size [l4|. However, they caution that 
the particular problems they studied have a high degree 
of structure that would also make them easy to solve 
classically. Nevertheless, the results show that the adi- 
abatic algorithm scales favorably at least on easy prob- 
lems. To show the gap is, at the very least, non- vanishing, 
Ruskai 0] provides a clever proof that the ground state 
of the instantaneous Hamiltonian must be unique. How- 
ever, as she points out, this tells us nothing about the 
magnitude of the gap, and how it scales with problem 
size. A less encouraging result was obtained by van Dam, 
Mosca, and Vazirani who were able to construct a 
family of minimization problems for which they could 
prove an exponential lower bound on the running time of 
the (original) adiabatic algorithm on these problems. A 
subsequent paper by Farhi et al. challenged the in- 
evitability of such results by arguing that they might be 
circumvented by choosing a different interpolation path 
between the initial Hamiltonian and the one encoding 
the problem to be solved. To date, the most sophisti- 
cated analysis of the running time of the adiabatic algo- 
rithm on NP-complete problems was provided by Roland 
and Cerf . They found that by nesting one quantum 
adiabatic search algorithm within another, one can solve 
NP-complctc problems more efficiently than naive use of 
an adiabatic version of Grover's algorithm. Nevertheless, 
the run-time scaling is still exponential in problem size, 
albeit better than what is possible classically. 



IV. APPLICABILITY OF RANDOM MATRIX 
THEORY 

Random Matrix Theory is a statistical approach to 
Hamiltonian systems that are otherwise analytically in- 
tractable. For example, RMT focuses on universal 
model-independent properties of the system under study, 
such as the distribution of the spectral fluctuations. 
Many superficially different physical systems are found 
to have distributions of spectral fluctuations that fall into 
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just a handful of categories. Such a characterization of 
spectra originated in the context of nuclear physics [l9l | 
and was applied later to complex many body systems and 
quantum systems having a chaotic classical analog pcj . 
Once the distribution of spectral fluctuations of a phys- 
ical system has been identified (and deemed to be irreg- 
ular), RMT can be applied to make predictions about 
properties of interest, such as transition rates between 
different levels. For example, the problem of estimating 
transition rates has been examined in the context of nu- 
clear dissipation, and the use of the LZ transition as a 
mechanism for nuclear dissipation was suggested origi- 
nally by Hill and Wheeler [21|. The combination of the 
LZ transition probability with the RMT statistical ap- 
proach was examined by Wilkinson j22, H3, whose re- 
sults we apply to the current problem. Although the 
LZ assumption is reasonable for adiabatic systems |25j |. 
the mechanism for dissipation in complex spectra con- 
tinues to be investigated in the context of RMT, with 
more recent approaches using a non-LZ, propagator ap- 
proach [26j . 

To determine the applicability of RMT, the regular- 
ity of the entire spectrum at each adiabatic interpolation 
point is measured by the Brody parameter (defined be- 
low). The nearest- neighbor spacing distribution is the 
most common measure of spectral regularity in quantum 
systems. This measure is quantified by the Brody pa- 
rameter q that interpolates between a regular Poisson 
spectrum (q = 0) and an irregular (quantum chaotic) 
Wigner distribution (q = 1) |10| . A renormalized spec- 
trum with a Brody parameter q is characterized by the 
following NNS probability distribution for spacing level 
spacing 5: 

/ o _i_ \ "1 1 +9 
p q (S) = (1 + q)(35« exp(-/3<S 1+ «) , (3 = T ^ 

(6) 

The form of this distribution, for different values of the 
Brody parameter, is shown in Fig. [21 

Irregular RMT spectra are characterized by an abun- 
dance of avoided level crossings, and a lack of level degen- 
eracies. We now proceed to determine whether or not the 
distribution of spectral fluctuations is anywhere irregular 
during the interpolation process. 



V. SPECTRAL FLUCTUATION EXPERIMENTS 

Our first task is to determine an appropriate ensem- 
ble of random soluble 3-SAT problem instances to use. 
We would like to use computationally hard problem in- 
stances, because we are most interested in assessing the 
scaling of the failure rate of the adiabatic algorithm on 
hard problems. Hard, in this sense, is a relative term. 
When solving random instances of soluble 3-SAT prob- 
lems having n variables and m clauses, typically the hard- 
est instances are encountered at a critical value of the 
clause to variable ratio m/n. For the 3-SAT problem as 
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FIG. 2: Nearest neighbor probability distribution p q (5) for 
eigenvalue spacings 5, in unfolded spectra for values of the 
Brody parameter q £ [0, 1]. When q — 0, the distribution re- 
sembles an exponential distribution, but changes to unimodal 
as q is increased. 



n — ► 00, the hardest instances are clustered around the 
ratio m/n w 4.2. However, small problem instances (hav- 
ing, say, n < 40), typically have a somewhat displaced 
transition point. Fig. 01 shows the mean computational 
cost of solving 3-SAT problems containing from n = 8 
to 11 = 50 variables using either the Davis-Putnam (DP) 
algorithm Q or the GSAT algorithm 0. Regardless of 
the algorithm used, an easy-hard-easy pattern is appar- 
ent when the number of clauses is increased at fixed num- 
ber of variables. In the limit of infinite problem size, the 
easy and hard instances are separated by a phase tran- 
sition (see, e.g., [2!|). But, as can be seen from Fig. [21 
the location of the phase transition point is extremely ill- 
defined for problem instances having n < 20. Thus infer- 
ring any reliable cost-scaling by extrapolating costs from 
such small instances would be exceedingly unreliable. Yet 
this is exactly what has been done in assessing the scal- 
ing of the quantum adiabatic algorithm from numerical 
experiments @|. As the simulation of the adiabatic al- 
gorithm solving an n variable 3-SAT problem involves 
2>i x 2 n -dimcnsional matrices, we cannot simulate very 
large cases. Hence, we begin by first determining the ac- 
tual location of the hardest problems for 3-SAT problems 
involving a more tractable n = 8 variables, rather than 
relying on the asymptotically known result whose appli- 
cability is suspect at small n. Specifically, we generated 
72,000 random 3-SAT problem instances all having n = 8 
variables, but with the number of clauses ranging from 
m = 8 to 80, corresponding to 1 < m/n < 10, and solved 
them using the GSAT algorithm. Each data point was 
computed from an average of 1000 problem instances. 
The results are shown in the lower portion of Fig. [3J Al- 



5 



s=0.06q=0. 



s = 0.22 q = . 



4000 
H 3000 

H 

nj 
u 

Q 2000 
1000 



var 




2 3 4 5 6 7 8 
Ratio of Clauses-to-Variables 



H 




SA 


80 




60 




40 


a 


te 


20 







=#= 


G 


-20 


n! 


Me 


-40 



2 3 4 5 6 7 8 
Ratio of Clauses-to-Variables 



3 . 5i 
2.5 

4 

0.5 



5 10 15 20 25 30 
s = 0.4 q = . 



L 



L 



5 10 15 20 25 30 
s=0.56q=0. 



L 



. ! 



5 10 15 20 25 30 
s=0.74q=0. 



5 10 15 20 25 30 
s = . 9 q = . 



1.5| 

1.25| 

. 75 
0.5 
0.25 




5 10 15 20 25 30 



FIG. 4: Distribution of eigenvalue gaps in the unfolded spec- 
tra of the instantaneous Hamiltonians during the course of 
the adiabatic algorithm solving easy instances of 3-SAT hav- 
ing n = 8 variables, and m = 4 clauses, as the interpolation 
parameter varies from to 1. As s increases, the Brody pa- 
rameter remains zero, implying that RMT is not applicable 
to easy problem instances. The horizontal axes are nearest 
neighbor eigenvalue spacings (NNS), expressed in tenths of 
the mean local spacing. 



FIG. 3: Upper plot shows the mean computational cost of 
solving n = 20, 40, and 50 variable 3-SAT problems using 
the Davis-Putnam algorithm (the data on the scaling of this 
algorithm is due to Bart Selman, see [2S|). The data show that 
for small problems (i.e., n < 50) the phase transition region 
is smeared out. Lower plot shows the mean computational 
cost of solving 1000 random instances of guaranteed soluble 
3-SAT problems having n = 8 variables, and m = 8 to 80 
clauses, using the GSAT algorithm. Again, for such small- 
sized problems the region of hard problem instances (relative 
to other instances) is quite spread out, but roughly centered 
on m/n — 6, rather than the asymptotic value of 4.2. 



though the difference between easy cases and hard cases 
is not as pronounced as it is for much larger n, neverthe- 
less, the data suggests that problem instances centered 
around m/n = 6 will yield relatively hard cases [13|. We 
will use such problem instances to create the ensemble 
we need in our numerical studies of the distribution of 
spectral fluctuations of the instantaneous Hamiltonians 
encountered during the interpolation phase of the adia- 
batic algorithm. 

Next we turn our attention to the global spectral prop- 
erties of the instantaneous Hamiltonians encountered in 
the quantum adiabatic algorithm for easy and hard prob- 
lem instances. Specifically, we obtain the NNS distri- 



bution of the adiabatic Hamiltonian H(s) after renor- 
malization of the spectrum to unit average local level 
density 0, @ for an ensemble of easy cases, and for an 
ensemble of hard cases. The regularity of the spectrum 
is determined at each point in the adiabatic evolution 
by fitting it to ©, and obtaining the Brody parameter 
q(s) of the NNS distribution at that point. We begin 
by determining the spectral distribution of an ensemble 
of easy problems as the interpolation parameter in the 
quantum adiabatic algorithm ranges from s = to s = 
1, for instances of soluble 3-SAT having n = 8 variables 
and m = 4 clauses, i.e., problems for which m/n = 0.5. 
Each histogram in Fig. 0] is based on the spectral behav- 
ior of an ensemble of 20 problem instances having fixed 
values of n and m. For the easy problems, the spectrum 
of each instantaneous Hamiltonian conforms to a Pois- 
son (regular) spectral fluctuation distribution, and small 
Brody parameter, q = 0. Hence, we conclude that RMT 
is not applicable in the easy region, i.e., for m/n < 5 
(n=8). 

In contrast, Fig. [3] shows the eigenvalue gap fluctua- 
tions of instantaneous Hamiltonians induced from "hard" 
instances of soluble 3-SAT with n = 8 variables and 
m = 48 clauses. Again, our ensemble averages over 20 in- 
stances of a fixed n and m. For m/n = 6 (hard problems) 
we observe Poisson behavior for s < 0.5, but for s > 0.5 
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FIG. 5: Distribution of eigenvalue gaps in the unfolded spec- 
tra of the instantaneous Hamiltonians during the course of the 
adiabatic algorithm solving hard instances of 3-SAT having n 
= 8 variables, and m = 48 clauses, as the interpolation param- 
eter varies from to 1. As s increases, the Brody parameter 
varies from to 0.49, and back to 0. The horizontal axes 
are nearest neighbor eigenvalue spacings (NNS), expressed in 
tenths of the mean local spacing. 



the spectra become increasingly irregular and the Brody 
parameter becomes significant. In other words, the in- 
stantaneous Hamiltonians induced by random, hard 3- 
SAT instances, appear to have a qualitatively different 
spectrum from those of easy problems of the same size. In 
particular, at a certain point in the interpolation process 
between the initial and final Hamiltonians, the spectrum 
becomes irregular, and the NNS distribution resembles a 
Wigner distribution with a relatively large Brody param- 
eter. Here, RMT can be applied to estimate transition 
rates between levels. 

Finally, one can repeat these experiments for problems 
having n = 8 variables at a clause to variable ratio of 
m/n = 9 (easy region again - data not shown). Here, we 
observe results similar to m/n = 6, but with a slightly 
smaller maximum Brody parameter. We attribute this 
to the fact that m/n = 9 is easier to solve than m/n = 6, 
although for n — 8 variables, this difference is slight. 

Fig. E| summarizes the Brody parameter as a function 
of the interpolation parameter s for easy and hard prob- 
lems with n = 8 variables, and from 4 to 72 clauses. The 
critical question in deciding if the quantum adiabatic al- 
gorithm can be completed in polynomial time is whether 
such a (fast) interpolation would induce level transitions. 
If they do occur, the system will not reside in the ground 
state of Hi upon completion of the adiabatic path, and 
the algorithm will have failed to find the solution. In the 
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FIG. 6: Brody parameter as a function of the interpolation 
parameter s,for easy and hard instances of 3-SAT with n = 8 
variables from m — 4 to 72 clauses. Note that the largest val- 
ues of the Brody parameter coincides with the hardest prob- 
lems. For soluble 3-SAT problems of size n = 8, these are 
found around m/n = 6 rather than at m/n = 4.2, which is 
the asymptotically valid transition point for 3-SAT as n — > oo. 



next section, we compute the probability that level tran- 
sitions occur during the course of the interpolation from 
Ho to Hi using RMT, for problems of a given degree of 
difficulty. 



VI. RANDOM MATRIX ANALYSIS OF 
QUANTUM ADIABATIC ALGORITHM 

Previous analyses of the adiabatic algorithm (summa- 
rized in Section placed the greatest significance on 
the scaling of the E± — Eq gap with increasing problem 
size. However, as Equation (JHJ shows, it is the ratio of 
the matrix element (dH/dt)i^ to the square of the mini- 
mum gap that determines whether the adiabatic theorem 
applies. Here, we calculate instead the probability of a 
transition from the ground state, which is a reliable proxy 
for the failure rate of the adiabatic algorithm. In regions 
where the Brody parameter is significant, we assume that 
any non-adiabatic transitions arc of Landau-Zcncr type, 
i.e., confined to adjacent levels at avoided crossings where 
the energy levels assume the geometry of convergent hy- 
perbolae [ll|] (see Fig.nj. When approached in this man- 
ner, the probability for a single transition, anywhere in 
the spectrum, can be parameterized by the minimum gap 
size AE, the difference in the asymptotic slopes Am, and 
the rate of change of the adiabatic evolution parameter 
ds/dt = s. Specifically, the transition probability is [Tlj 



-2-7T7 
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1 AE 
Ah \Am\i 



(7) 



Typical values for the parameters AE and Am will vary 
with the difficulty of the problem instance being solved 
(reflected by the clause to variable ratio m/n) as well as 
with the interpolation parameter. 

In order to exploit this transition probability to predict 
the rate of transition from the ground state, we need to 
verify that the Eq — E\ gap fluctuations follow the same 
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distribution as those in the body of the spectrum. In 
Appendix A, we show that the distribution of gap fluctu- 
ations is characterized by a Brody parameter comparable 
to the typical fluctuations in the body of the spectrum, 
establishing this point. To incorporate the model de- 
pendence (non-universality) of the problem, the average 
local level density of the spectral region of interest (here 
the levels) is included. Thus, the RMT/LZ calculation 
results in an ensemble-averaged transition rate with the 
LZ process as the principle transition mechanism. 

The transition rate from the ith eigenstate of the in- 
stantaneous Hamiltonian with energy Ei, dPi/dt, can be 
written in terms of the second moment of the occupation 
probability distribution Pi{t) because it essentially is 
the rate of diffusion of the occupation probability. As dis- 
cussed in the introduction, the appropriate RMT ensem- 
ble to use for physical systems with exponential densities 
of states is the embedded GOE, consisting of real, orthog- 
onal matrices having Gaussian-distributed random ma- 
trix elements with at most two-body interactions. Since 
in the calculation of transition rates in RMT only the 
fluctuation properties of the spectrum (rather than the 
density of states proper) enter (2^, we can safely substi- 
tute a GOE to obtain the probability to transition from 
the ground state 



dP 
~dt 



oc^/yi^i 3 / 2 



(8) 



where p is the level density of Eq — E\ levels averaged 
over an ensemble of problem instances having the same 
clause to variable ratio, and a is the typical size of the 
asymptotic slopes of the LZ avoided level crossings in 
the region of interest. Indeed, as shown in Ref. [24||, em- 
bedded GOEs give rise to NNS distributions very sim- 
ilar to those arising in a GOE, except that the maxi- 
mal Brody parameter is limited to q w 0.8. Note that 
this is also implied by a numerical calculation p5| of the 
dissipation rate of probability distributions of physical 
instantaneous Hamiltonians (i.e., those drawn from and 
embedded GOE), which agree with the scaling in Eq. ||HJ). 

We emphasize that Eq. [S] expresses the transition rate 
from the ground state in terms of the average Eq — E\ 
level density, rather than the minimum gap (i.e., the max- 
imum level density). The average level density is the 
more relevant parameter for assessing the typical behav- 
ior of the adiabatic algorithm. As the adiabatic algorithm 
fails if the system transitions from the ground state, Eq.|S| 
can be interpreted as the average failure rate of the adia- 
batic algorithm when driven with uniform interpolation 
velocity ds/dt. 

Finally, we need to determine the scaling of the transi- 
tion rate from the ground state (i.e., the failure rate) with 
increasing problem size. We do this in two stages: first we 
show how the transition rate must scale with p in order 
to keep the transition rate bounded, then we show how 
p scales with problem size. Although the values of p and 
<7 are model-dependent quantities, we note that a is the 
characteristic slope in the energy-parameter space (i.e., 



a ks dE/ds). To make explicit the dependence on level 
density, we work with unfolded energies {Ei = A x e.;) 
where A is the mean level spacing (A = 1/p) and are 
the unfolded energy levels having mean level spacing one. 
Under this transformation, er 3 / 2 — > (A x <r) 3 / 2 , and we 
therefore write 



dP 

HE 



=oca 3 /y/ 2 |s| 3 / 2 



(9) 



In general, while the unfolded level slope a and the av- 
erage level density p are model-dependent quantities [27J , 
Eq. © nevertheless exhibits the explicit dependence of 
the transition rate on p. We can establish a lower bound 
on the interpolation time required to evolve the system 
through the irregular region by noting first that 



dP . dP,. 
-^( S ,p)>^(s,p mia ) 



(10) 



where we defined p m i n = min s p(s) and where the min- 
imization is carried out over only those values of s that 
give rise to an irregular spectrum. It is important to note 
that p m in is not the minimum of a particular problem in- 
stance; it is an average level density at a particular value 
of s, whose value is constant for a given problem parame- 
ter set. The failure probability during an evolution path 
of length T is then bounded by 



P > const5- 3/2 As 3/2 



T 



1/2 



(11) 



To ensure a given transition probability over a given 
range As in an irregular region, the interpolation time T 
must scale as p m ia- If we can now estimate how p mm (an 
average quantity) scales with increasing problem size, we 
can estimate how the time needed to complete a partic- 
ular part of the adiabatic algorithm must scale in order 
to keep the transition rate from the ground state small. 
This part is precisely the region where one would need to 
go most slowly to avoid an unwanted transition from the 
ground state. What can we say about the scaling of p m i n 
with problem size? As pointed out by Ruskai [15| . in re- 
gions characterized by a lack of level degeneracies (such 
as irregular spectral regions) the interpolating Hamilto- 
nian H (s) must fit 2™ eigenvalues into a range that is 
polynomial in n, and consequently the level density must 
scale exponentially with problem size. Our numerical 
simulations provide strong evidence for the existence of 
such irregular spectral regions by the fact that we obtain 
significant values of the Brody parameter for distribu- 
tions averaged over instances that have a sizable degree 
of difficulty (as measured by their clause to variable ratio. 
In other words, difficult problem instances display spec- 
tra characterized by a large Brody parameter and few 
level degeneracies, and therefore show a lack of symme- 
try Thus, with the final assumption that non-degencrate 
spectra (irregular spectra having a large Brody param- 
eter) have a minimum average level density that scales 
exponentially with problem size, we find the transition 
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rate out of the ground state (the algorithmic failure 
rate), and thus that the adiabatic evolution time cannot 
scale polynomially with problem size. 

The information regarding the parameter regions hav- 
ing significant transition probability may be used in a 
variable-speed approach to the adiabatic quantum algo- 
rithm in a manner similar to fl8l | . In this way, a speedup 
over classical algorithms for solving NP-complete prob- 
lems still appears possible. We wish to emphasize that 
an d the entire RMT analysis, holds only in those 
regions of the interpolation process where we find irreg- 
ular spectra, i.e., non-trivial Brody parameters. In par- 
ticular, the formula does not apply to regions of the in- 
terpolation having a small Brody parameter, arising typ- 
ically in the very early and very late stages of a hard 
problem, or at all stages of an easy problem. There- 
fore, our formula is not inconsistent with previous results 
(e.g., in Section HTT1 showing that the adiabatic algorithm 
scales favorably on easy problems) , because in such cases 
Eq. (fTT|) is invalid. 



APPENDIX 



Figure 7 shows the ground state level distribution of 
Eq — Ei gaps for 1000 random problem instances in the 
hard region for problems of size n = 8. The horizon- 
tal axis is the eigenvalue spacing in tenths of the mean 
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VII. CONCLUSION 

Our approach to the analysis of Adiabatic Quantum 
Computation is statistical in nature, and relies on the 
applicability of Random Matrix Theory to large Hamil- 
tonians. With such an approach, we examine the be- 
havior of an averaged failure rate and interpolation time 
when averaged over an ensemble of problem instances of 
a given difficulty. We find a spectral transition from or- 
derly to disorderly as problem difficulty is increased. Our 
analysis of these results suggests using the Landau-Zener 
transition probability in a statistical RMT approach. For 
those regions where RMT applies, degeneracies are lack- 
ing and we show that average failure rate, and average 
interpolation time do not scale polynomially with prob- 
lem size. Nevertheless, the quantum adiabatic algorithm 
appears to be more efficient than any known classical al- 
gorithm for solving NP-complete problems, with a speed- 
up commensurate with Grover's search algorithm. More- 
over, the insights we gain from the spectral fluctuation 
analysis may be used to design a variable interpolation 
rate, which has the potential to improve the adiabatic 
algorithm over its performance using a uniform interpo- 
lation rate. 

Finally, we believe the connection shown here between 
difficult computational problems and physical spectral 
irregularity is a powerful path towards understanding and 
classifying physical approaches to computation. 



FIG. 7: NNS distribution of E - £i gaps for 1000 random 
problem instances in the hard region for problems of size n = 
8. The horizontal axis is the eigenvalue spacing in tenths of 
the mean level spacing. 



level spacing. The distribution shows that spectral ir- 
regularity definitely extends to the lowest lying levels for 
hard instances of 3-SAT problems, implying that RMT 
is applicable to hard instances of 3-SAT across the en- 
tire spectrum of levels, including the Eq — E\ gap, which 
is the relevant gap for assessing the scaling of the quan- 
tum adiabatic algorithm. Therefore, our use of RMT to 
characterize the transition rate out of the ground state, 
and hence to estimate the cost scaling of the quantum 
adiabatic algorithm, is valid. 
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